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We present the non-forward quark matrix elements of operators with one and two covariant derivatives needed 
for the renormalisation of the first and second moments of generalised parton distributions in one-loop lattice 
perturbation theory using clover fermions. For some representations of the hypercubic group commonly used 
in simulations we define the sets of possible mixing operators and compute the one-loop mixing matrices of 
renormalisation factors. Tadpole improvement is applied to the results and some numerical examples are presented. 
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1. INTRODUCTION 

Generalised parton distributions (GPDs) have 
become a focus of both experimental and theo- 
retical studies in hadron physics (for an exten- 
sive up-to-date review see P). They allow a 
paramctrisation of a large class of hadronic cor- 
relators, including e.g. form factors and the ordi- 
nary parton distribution functions. Thus GPDs 
provide a solid formal basis to connect informa- 
tion from various inclusive, semi-inclusive and 
exclusive reactions in an efficient, unambiguous 
manner. Furthermore they give access to physical 
quantities which cannot be directly determined 
in experiments, like e.g. the orbital angular mo- 
mentum of quarks and gluons in a nucleon (for a 
chosen specific scheme) and the spatial distribu- 
tion of the energy or spin density of a fast moving 



*Talk given by A. Schiller at the Workshop on Computa- 
tional Hadron Physics, Nicosia, September 2005. 



hadron in the transverse plane. Since the struc- 
ture of GPDs is rather complicated a direct ex- 
perimental access is limited. Therefore, comple- 
mentary information channels have to be opened 
up. One major source is lattice QCD |2I3I4I5Io] . 

Recently [7] we have calculated the non- 
forward matrix elements in one-loop lattice per- 
turbation theory for the Wilson fermion action 
needed for the renormalisation of the second mo- 
ments of GPDs. From these results the renormal- 
isation factors for various representations of the 
hypercubic group have been derived. 

Here we present some new results for oper- 
ators with two covariant derivatives using the 
Shcikholcslami-Wohlcrt (clover) action [H], which 
leads to 0(a) improved quark-quark-gluon and 
quark-quark-gluon-gluon vertices in the Fcynman 
rules. Since in current numerical simulations the 
operators for the second moment of GPDs are 
not improved, we ignore such a possible addi- 
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tional improvement. Note that the O(a) improve- 
ment for operators with one covariant derivative 
is known [§]. 

We consider the Wilson gauge action and clover 
fermions with the fermionic action Ssw,f [HI (f° r 
dimcnsionful masslcss fcrmion fields tJj(x)) 

S S W,F = W^(x)V(x) (1) 

x 

a 3 

+ tp(x + ap)(r + ln)Ul ifJt il)(x)] 

Here a denotes the lattice spacing and the sums 
run over all lattice sites x and directions /i, v 
(all other indices are suppressed). F^° VCT is the 
standard "clover-leaf" form of the lattice field 
strength and a^ v = i/2[<y fl ,j„}. 

In a perturbative calculation the operators to 
be investigated are sandwiched between off-shell 
quark states with 4- momenta p and p' . Our calcu- 
lations are performed in Feynman gauge, the final 
numbers are presented for the Wilson parameter 
r = 1, leaving the value of c sw free. 

One-link quark operators with clover fermions 
have been discussed in [§] for forward matrix 
elements. The renormalisation constants found 
there for a given representation of the hypcrcubic 
group H(4) and charge conjugation parity C can 
be used in the non-forward case as well. Addi- 
tionally, in the case of GPDs "transversity" op- 
erators have to be taken into account. We will 
collect here all results for completeness. 

It is well known that operators with two or 
more covariant derivatives may mix under H(4): 
the one- and higher-loop structures differ in gen- 
eral from that of the Born term and multi- 
plicative renormalisation may get lost. In addi- 
tion, for non-forward matrix elements also opera- 
tors with ordinary (external) derivatives can con- 
tribute making the mixing problem more compli- 
cated. To find the possible candidates for mixing 
one has to define those operators which belong to 
the same irreducible representation under H(4) 
and have the same charge conjugation parity. 



We define renormalised operators Of by 
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Of{n)=Y,Zf h {a,n)O k {a) 



(2) 



k=l 



where S denotes the renormalisation scheme and 
N is the number of operators which mix in one- 
loop. Zf k (a,fi) are the renormalisation constants 
connecting the lattice operator Ok (a) with the 
renormalised operator Of {fx) at scale fi. We 
present the renormalisation constants in the MS 
scheme following 0. 

2. OERATORS AND MIXING 

We consider operators with up to two covariant 

symmetric lattice derivatives D=D — D and ex- 
ternal ordinary derivatives d needed for the cho- 
sen representations of interest for the first and 
second moment of GPDs. The standard reali- 
sation of the covariant derivatives acting to the 
right and to the left is used: 

D» ift(aO = ±- x (3) 



2a 



U x ^ tj)(x + ap) - Ul_ a ~ tlx 4'{x - ap) 
tp(x) D„= — x 



2a 

ip(x + ap) Ul ^ - tp(x - ap) U x -ap.,v, 



The external ordinary derivative is taken as 

% • • • V>) (x) = - x 
(ip - ■■ if) (x + ap) - (ip ■ ■ • if) (x) 



(4) 



(5) 



The number of derivatives appearing in the op- 
erators is indicated by superscripts D and d, re- 
spectively. Quark operators with one derivative 
are given by 
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5,D 
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T,D 

flUUJ 



o 



T,d 



"07m 75 D v ^ , 



(6) 
(7) 
(8) 
(9) 
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The operator © is a transversity operator anti- 
symmetric in its first two indices which is of inter- 
est for GPDs, operators JHJl and © contribute as 
lower dimensional operators to mixing in certain 
representations of the second moment of GPDs. 

As operators with two derivatives we consider 
here 



O 



DD 

flUUJ 



1 



T^lu DvDu 1p ■ 



O 



dD 

flUUJ 



o 



dd 

flUUJ 



(10) 



In addition, spin-dependent and "transversity" 
operators have to be considered when discussing 
all possible representations. They are roughly ob- 
tained by replacing 7^ by 7^75 and a^r, respec- 
tively. 

To define the various representations with given 
C we use the following short-hand notations 



0{ Vl V 2 V 3 } 
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1 

6 
O 



(0-V\V<2 
V\ VI V 3 
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\V1V2V 3 \\ 



O - o 
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(<!/! V 2 V 3 )) 



(0 
2 O 

<s V 2 V\ V 3 ; 

a, + a 
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Let us denote an irreducible representation of the 
hypercubic group H(4) by rjp with dimension Z 
(k labels incquivalent representations of the same 
dimension) and a given charge conjugation parity 
C by ±1. 

For the first moments we choose the following 
representations presented in Tabled (for the nota- 
tion and a detailed discussion of the transforma- 
tion under H(4) see ^Hl)- They are renormaliscd 
multiplicativcly. 

For the second moments we consider in this 
contribution the following mixing cases ( the de- 
tails and additional operators will be presented 
elsewhere 



Table 1 

Operators and their transformation under the hy- 
percubic group. 



Operator 



n 



— 



C 



u {n} 




+ 1 


og - 1 (oft + og + og) 




+ 1 


U {14} 


r (6) 


-1 


otf-i(oif + o 5 2 f + otf) 


r O) 
T 4 


-1 


T,D 
U «124)> 




+ 1 


Q T,D _ T,D 
U ({122)) U «133» 


r« 


+ 1 


n T ' D -O t ' d 

W \\122\\ || 133 1| 


r<* 


-1 



Representation Tjj , C = — 1 with operators 



U {124,} > ^{124} 



(11) 



(8) 

Representation t\ , C = — 1 with 

O 



Ox 

o 3 
o 7 



1 



U {114} 2 \ U \2 



DD 
{224} 



U {114} " 

n DD 
U «H4» 



I ( 



'{224} 







2 !, U «224» 



{334} J i 
{334} J i 
U <(334)) 



U «114» - 2 (, U «224)> 



°f<334»j 



5,9D 
1 1 2 13 1 1 ' 
5.DD 
||213|| ' 



o 6 = o 



«213» 



(12) 



,9 _ 1 / n T,d 
11 2 V 422 



^411 



o. 



433 



3. ONE-LOOP CALCULATION 

We calculate the non-forward matrix elements 
of the operators in one-loop lattice perturba- 
tion theory in the infinite volume limit following 
Kawai et al. [T^J. Details of the computational 
procedure are given in |7j. 

In lattice momentum space operators with non- 
zero momentum transfer q are realised by apply- 
ing the lattice momentum transfer q at the lattice 
position x or at the "position centre" x + e.g. 
for an operator with one covariant derivative we 
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have the two possibilities 

(*M( ff ) = j££{' 

x v 

[tjj(x)U x> ^(x + afi) - i>(x + afiU* i/>(x)] . (13) 



e ig-x _|_ gig-(a:+o/l) 
2Qiq-(x+afi/2) 



Eq. I|13|) basically defines the Feynman rules for 
the operators in lattice perturbation theory. As 
an example we get for the operator to order 

0(g°) (g is the bare gauge coupling): 



O 



°Z{p',p) = i>(p')iMp) x 

1 . a(p + p') v 1 . a(p + p') u 

- sin sin x 

a 2 a 2 



(14) 



cos 



a(p-p')- 



COS 



In the following sections we denote the up- 
per/lower realisations by supercripts I /II. 

The contributing one-loop diagrams for the self 
energy and the amputated Green (vertex) func- 
tions are shown in Figures ^ an d 121 (filled black 
circles indicate the place of the operator inser- 
tions): 






Figure 2. Amputated Green function diagrams. 



Table 2 

Finite contributions B(c sw ) for the first moments. 



Representation 



B(c sw ) 



4 6) 


1.280- 


3.873 c sw - 


0.678 cL 




2.562 - 


3.970 Csw 


1.040 cL 


r (6) 
r 4 


0.345 - 


1.359 c sw - 


1.893 cL 


r O) 
r 4 


0.167- 


1.249 c sw - 


1.998 cL 




13.169 4 


- 2.675 c sw - 


1.494 cL 


rf (C = ±1) 


12.804 -t 


- 2.624 c sw - 


1.430 cL 



Figure 1. Quark self energy diagrams. 



3.1. First Moment 

Since mixing is absent, we omit the matrix no- 
tation for the renormalisation constants and use 
the general form {Cf = 4/3, gn is the renor- 
malised coupling) 



Z(a/i) 



1 



9lC F 
16tt 2 



[ 7 ln(aV)+£(c™)] (15) 



with the anomalous dimensions 7 = 8/3 for the 
first four operators and 7 = 2 for the last three 
operators in Table For the operators of that 
Table we get the finite contributions shown in Ta- 
ble □ 



3.2. Second Moment 

We present the matrix of renormalisation con- 
stants in the generic form 



with 



9rCf 
16tt 2 



7<iMaV)+-B{j n) (c™) (16) 



Bif(c sw )=B, 



(0,m) _,_ g(l,m) 



1 J 



„ B (2,m) 2 



The superscript (m) with m = I, II distinguishes 
the realisations I and II of the covariant deriva- 
tives i(i"3l) . 

Representation C = — 1 

For this representation the operators Hlfl mix. 
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The corresponding 2 x 2-mixing matrices are 



7jfc 



2f> 
6 





H jk ( c ™>) 



(17) 



-11.563 0.024 
20.618 





0.984 




-0.255 
4.746 

0.016 
0.543 



.(18) 



the lower dimensional operator Os given in (|12() . 
The perturbative mixing result is 

n I — QlfiF 1 , tree 

^ll/a-part" lg7r 2 a U S 



(20) 



(-0.518 + 0.0832 c s 



0.00983 cjj, 



but a nonperturbative subtraction from the ma- 
trix element of 0\ is required to obtain reliable 
numbers. 

4. TADPOLE IMPROVEMENT AND 
SOME NUMERICAL EXAMPLES 



In the matrix -B^ the mixing between the op- 
erators OP^ 4 x and 0{f 2 4} ^ s verv small. Thus it 
may be justified to neglect the mixing in practical 
applications. 



Representation 



C= -1 



We consider the mixing (|12f> of the operators 
having the same dimension first. These are the 
operators Oi . . .Oj in (p~5f) . To one-loop accuracy 
the operator O7 does not contribute and we have 
to consider the following mixing set: 

{0!, 02, 03, 04, 05, 6 }. 



The anomalous dimension matrix is 



( 25 
' 6 






V 







5 

6 



I 









\ 


2 


-2 

2 
3 



(19) 



and the finite parts of the mixing matrix are given 
in TableEl (in cases of doublets the upper number 
belongs to type I, the lower to type II of realisa- 
tion of lattice covariant derivative) . 

Using lattice perturbation theory to one-loop, 
l/a terms may appear when calculating the ma- 
trix elements of the operators with two covariant 
derivatives. Such terms are potentially danger- 
ous because of the power-law divergence in the 
continuum limit. Considering the representation 
Tj 4 ', a potential mixing is absent. On the con- 
trary, we get mixing for operator 0\ of r\ with 



Since many results of (naive) lattice pertur- 
bation theory are in bad agreement with their 
numerical counterparts, it has has been pro- 
posed to rearrange the (naive) lattice pertur- 
bative series. This rearrangement is performed 
using the variable uq (the mean field value of the 
link), e.g. defined from the measured value of the 
plaquette at a given coupling 

u = (^TvU n )i. (21) 

In case of mixing the tadpole improvement pro- 
cedure proceeds as follows. By scaling the link 
variables Up with uq 



U^x) = Uq 



Uq 



U U,J,{X) 



(22) 



the amputated Green function for operator 
with n covariant derivatives A^ 1 takes the form 



A 



(n) 



u%A%\u„(x)) 



(23) 



A@ (U^x)) is expected to have a better converg- 
ing perturbative expansion. Up to order g 2 we 
obtain for the Wilson gauge action, labelling the 
operators by i and the corresponding number of 
covariant derivatives by n,, 



^ ] (U,(x)) 



A<%L(U»(x)) 



pert 

9 2 C F 

16 7T 2 



n t ir 2 + 0(g 4 ) 



A 



(rii ,tree) 



9 2 C, 

16 TP 

-0( 5 4 )) 



^2 w lk k 



(nfc ,tree) 
k 



k=l 



(24) 



6 



Tabic 3 



(8) 

Finite parts of the mixing matrix for operators of representation r x , C 

( 



R (o,/,//) 



D 



jk 



D 



-12.127 


3.306 





2.922 


0.333 





-0.982 


0.371 







1.491 
-2.737 

20.618 

-8.015 
18.184 



3.264 
3.264 

-0.213 
-0.686 

4.746 

-0.766 
-0.055 



-1.441 
-1.441 

-0.078 
-0.101 

-0.543 

-0.551 
0.215 



1.416 
1.416 



0.368 


-14.852 







0.015 
0.994 







4.431 
-4.302 

20.618 





-0.033 


2.152 





-0.029 


-1.707 







0.015 
0.173 







1.206 
0.970 

4.746 



0.035 
0.042 



0.371 
0.116 

-0.543 





-1. 

0.016 0.150 



-0.928 0.738 



0.350 0.015 

0.005 0.360 / 

\ 



-0.019 0.057 


-1.758 2.298 



1.648 0.866 

0.289 2.225 / 



-0.001 




0.007 




-0.443 0.103 





-1.703 
0.189 





0.568 
-1.325 ) 



where the Wik denote the mixing weights. From 
(|24|l it becomes clear that in one-loop only the di- 
agonal terms in the mixing matrix get a shift pro- 
portional to n.;7r 2 . An external ordinary deriva- 
tive (d) does not provide a factor of u . Tak- 
ing into account the mean field value for the 
wave function renormalisation constant for mass- 
less Wilson fcrmions Zfjf^ ilson = u we get the 
tadpole improved matrix of renormalisation con- 
stants in the form 



(l - f^§(^ - 1) tt 2 h + 0(.g 4 )) Za . (25) 



Additionally, one has to replace the parameters g 
and c sw by their boosted counterparts 

9ti = 9 2 u o i > c " = c ™ u o ■ ( 26 ) 

Putting (US!), (HSJ and together we obtain 
for the tadpole improved renormalisation mixing 
matrix in one-loop order 

z Ti,(m) i_ ni / ffr/CV 

(7«InW)+^ /,(w,) (^))) (27) 

with 

<' (m) ( c ") = B$*{<%) + (»* - l > 2 ^ ■ (28) 



r 



Let us demonstrate the effect of tadpole im- 
provement by some numerical examples. We 
choose a = 1/fx, (3 — 6, u = 0.8778 and 
c S w = 1 + 0(g 2 ) For the first moments the 
only effect consists in replacing c sw by and 
9R by gn in (|15|l and in Table |5J For the repre- 
sentation ri 6 ^ we get 



Z = 1.02* 



Z TI = 1.023. 



(29) 



For the second moments we consider the simple 
mixing Of£ 4} <-> C^f 24} CU) first. Without tad- 
pole improvement we obtain the mixing matrix 



Zij — 



1.081 0.002 
0.790 



The tadpole improved result is 



zT/ = 



1.142 0.002 
0.707 



(30) 



(31) 



It might be instructive to compare the one- 
loop corrections for the renormalisation con- 
stants: b!-™\csw) for the unimproved case ifTTfl l 
and 5y (c^J) for the tadpole improved case 
(EZJ). We get 



Bi 



and 



B 



TI 



-9.649 




-0.209 




-0.247 
24.821 



-0.177 
12.035 



(32) 



(33) 



We observe that in agreement with the improve- 
ment aims the diagonal one-loop contributions 
are reduced. 

For the representation t[ 8 \ C = — 1 with the 
mixing of operators 0\ . . . Oq we obtain for the 
unimproved/improved mixing matrices (choosing 
m = I) the numbers given in Table 01 

5. SUMMARY 

Within the framework of lattice QCD with 
clover improved Wilson fermions and Wilson's 
plaquette action for the gauge fields we have cal- 
culated the one-loop quark matrix elements of op- 
erators needed for the first two moments of GPDs 
and meson distribution amplitudes. From these 



we have determined the matrices of renormalisa- 
tion and mixing coefficients in the MS'-scheme. 

For the first moments of GPDs we can use the 
results from the first moments of structure func- 
tions. The results for the second moments extend 
the numbers obtained with Wilson fermions [?]• 
The general conclusions concerning the mixing 
properties remain unchanged. All sets which con- 
sist of one operator with two covariant derivatives 
D and one operator with two external derivatives 
d show very small mixing. The set discussed here 
with seven potential candidates l|12|l shows a more 
significant mixing. 

Moreover, taking 0\ from (|12|l as the opera- 
tor to be measured in a numerical simulation a 
mixing with a lower dimensional operator 0% ap- 
pears. This requires a nonperturbative subtrac- 
tion for Wilson or clover fermions. 

Using overlap fermions, such a mixing with a 
dangerous lower dimensional operator must be 
absent, since the mixing operators are of differ- 
ent chirality. 
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